Quantum corrections to leptogenesis from the gradient expansion 
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Using the closed-time-path formahsm we quantify gradient corrections to the kinetic equations 
for leptogenesis, that are neglected in the standard Boltzmann approach. In particular we show 
that an additional CP-violating source term arises, which is non-zero even when all species are in 
local thermal equilibrium. In the early universe it is proportional to the expansion rate and would 
vanish for static equilibrium configurations, in accordance with the Sakharov conditions. We find 
that for thermal leptogenesis in a standard cosmological background the additional source term is 
small. However, it can become the dominant source in the limit of ultra-strong washout. 
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I. INTRODUCTION 

Today the observed universe almost entirely consists 
of matter, i.e. is baryonically asymmetric. An attractive 
explanation of the observed asymmetryis provided by the 
baryogenesis via leptogenesis scenario [l| . In this scenario 
the Standard Model is supplemented by heavy Majorana 
neutrinos. The CP- and lepton- number violating decay 
of the heavy neutrinos produces a net lepton asymmetry. 
The rapid expansion of the universe ensures that this 
asymmetry is not washed out by the inverse decay and 
scattering processes 01- Finally, anomalous electroweak 
processes convert the generated lepton asymmetry to the 
observed baryon asymmetry [3|, Ij] . 

The computation of the asymmetry in terms of the 
neutrino masses and mixing parameters requires a mi- 
croscopic description of the out-of-equilibrium decay pro- 
cess. The canonical approach is based on Boltzmann 
equations, furnished with decay and scattering rates com- 
puted from the vacuum S-matrix elements. The size of 
various corrections to this approximation has been in- 
vestigated by including e.g. deviations from kinetic equi- 
librium, quantum statistical terms, thermal masses and 
matrix elements [5l-ll0|. 

Beyond this it is important to check the validity of 
the semi-classical treatment embodied by the Boltzmann 
approach. This is particularly relevant for leptogenesis 
where loop effects and unstable particles are essential for 
the generation of the asymmetry. Within the standard 
Boltzmann (bottom-up) approach, these typically give 
rise to double counting problems as well as ambiguities 
related to the application of equilibrium quantum field 
theory for the analysis of out-of-equilibrium processes. 

These ambiguities can be resolved within the nonequi- 
librium closed-time-path or Schwinger-Kcldysh formal- 
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ism. In contrast to the bottom-up approach this may 
be viewed as a top-down approach. Here, the underly- 
ing microscopic description is based on the full quantum- 
mechanical evolution equation for the expectation value 
of the lepton current. Starting from the latter, quantum- 
corrected Boltzmann-like kinetic equations can be de- 
rived, which are inherently free of the double-counting 
problem, and include medium corrections to the CP- 
violating parameters in a consistent way |ll| - [l3| , see also 
[Tj, [la]. Furthermore, a possible influence of off-shell 
and memory effects has also been studied in this ap- 
proach [il-a^. 

The reduction of the full quantum equations of motion 
to Boltzmann-like equation requires a so-called gradient 
expansion in powers of space-time gradients dxi^ [20|. 
Physically, it amounts to an expansion in powers of the 
ratio of the microscopic time-scale tmic ~ ^/Mi and the 
macroscopic scales tmac ^ l/ri,l/i?. This means, it 
requires that the decay rates P; and the cosmic expansion 
rate H are much smaller than the corresponding right- 
handed neutrino masses Mi. 

The standard Boltzmann treatment of leptogenesis 
relies on the zeroth order in the gradient expansion, 
i.e. neglects all gradient terms. Although the condition 
tmic "^ tmac IS typically well fulfilled in thermal lepto- 
genesis, higher gradient terms could still be important 
for the calculation of the asymmetry, since the latter de- 
pends on the tiny difference between particle and anti- 
particle interactions and densities. At zeroth order t he 
deviation from equilibrium, which is crucial according to 
the third Sakharov condition, is described by the non- 
equilibrium distribution function of the heavy Majorana 
neutrinos: Af^.{p) = /^, (p) - f^-ip)- Consequently, the 
CP-violating source term is proportional to A/^ . . 

The gradient terms capture another potential source 
for a deviation from thermal equilibrium, which is due 
to the time-dependence of the effective temperature T{t) 
of the thermal bath of leptons, quarks, gauge- and Higgs 
bosons. Consequently, we expect that the gradient CP- 
violating source term is proportional to T = dT/dt. 
In contrast to the standard source, it remains non-zero 



in the limit A/^^ — >■ 0, and therefore could be impor- 
tant even though it is suppressed by the Hubble scale, 
T ~ —HT. Furthermore, it is conceivable that leptoge- 
nesis or baryogenesis occurs in the early Universe simul- 
taneously with other non-equilibrium phenomena such 
as (p)reheating or phase transitions in which case gradi- 
ent contributions could be strongly enhanced. Note that 
such gradient terms are crucial for electroweak baryoge- 
nesis, see e.g. |2l| - |27| . Finally, on the formal level, it is 
well-known that including first-order gradient terms en- 
sures the validity of exact conservation laws of the kinetic 
equations j28| . 

The aim of this work is to quantify the leading CP- 
violating source term proportional to T arising from gra- 
dient corrections to the kinetic equations. In section HIl 
we briefly review the Boltzmann and the closed-time- 
path approaches and set up our notation. In the same 
section we discuss the gradient expansion and the Boltz- 
mann limit. Next we analyze the new source term and 
discuss the results. 



(i) As we demonstrate in section lllli the additional 
source term has a qualitatively new structure and 
does not vanish even if all species are in local ther- 
mal equilibrium. Furthermore, we demonstrate 
that it becomes important for very heavy Majorana 
neutrinos and ultra-strong washout. 



(ii) In section llVl we argue that the gradient corrections 
also modify the well-known standard source term. 

Finally we summarize our results and conclude in sec- 
tion El 



II. GENERATION OF AN ASYMMETRY IN 
THE CTP APPROACH 

In this section, we first review the standard Boltz- 
mann approach, and then the closed-time-path (CTP) 
approach for the generation oi a B — L asymmetry. 
In order to illustrate the effect of the gradient correc- 
tions, we consider a toy model which has also been used 
in [ll|, [12: l29l . However, the derivation and the structure 
of the results are generic, and similar gradient corrections 
will also be present in phenomenological scenarios such 
as thermal leptogenesis. The Lagrangian of the model 
reads 

C = ^df'iP^d^i^^ - ]^Ml^,^H + d^'bd^b - m'bb 

-^m'~^_y^M-^m, ^^l,2. (1) 

It can be considered as toy model for a generic baryo- 
genesis scenario in which the asymmetry is produced by 
the out-of-equilibrium decay of some heavy species and 
where the CF-asymmetry in the decay is induced by the 
one-loop contributions depicted in Fig.[TJ For example, 
in thermal leptogenesis the real scalar fields ipi model the 




FIG. 1: Tree-level and one-loop-level diagrams of the decay 
process "0^ — >■ bb. The gradient contribution of the vertex 
diagram is suppressed by higher powers of the couplings and 
is neglected in the following. 



heavy Majorana (s)neutrinos, whereas the complex field 
b represents the (s)leptons. In GUT baryogenesis the 
real scalar fields model the heavy bosons and the com- 
plex field b the baryons. In the following, we shall simply 
refer to b as toy-baryons and to ipi as toy-neutrinos. 



A. Boltzmann approach 

The standard approach is based on generalized Boltz- 
mann equations for the distribution functions fa {X, p) of 
on-shell particle species a [30, Isil , 

p'^V^faip) - C, - Cr" [1 + faip)] - C'riaip) , (2) 

where Va is the covariant derivative, Ca are collision in- 
tegrals comprising gain and loss terms, and we suppress 
the space-time coordinate X for brevity. The latter take 
decays, inverse decays and scatterings into account, with 
rates inferred from the S'-matrix {in-out formalism). If 
one considers only the decay and inverse decay processes, 
see Fig.[Tl then for the toy-baryons 

Cr"(p) -/ dUldUl{27T)H{k -p- q)\M\l^,, 

X U,{k)[l + fb{q)] , (3a) 



C'rip) =JdUldUli2nr5{k -p-q)\M\ 
X [1 + Um]h{q), 



bb^%lji 



(3b) 



where dHp = d^p/[{2'!TY2Ep\. For the antibaryons /(, is 
replaced by f^. CPT invariance implies for the in-out ma- 
trix elements: \M\l^^,, = \M\l^^^^ = \\g^?{l + e\n 

and \M\l^_^-^j^ ^ \M\l,_^^ = \\g,\\l - e^'), where e"- 
is the usual CP-violating parameter. As is well-known, 
certain scattering contributions bb <->■ bb also have to be 
taken into account for consistency within the standard 
Boltzmann approach, see below. 

The total baryon density is given by the time compo- 
nent of the baryon current: 



nBit) = V-'JdVjoit,^) 
= rib- ni 



j (27r)3 yJb 



h] 



(4) 



To stress the analogy with phenomenological models, we 
consider the difference ub — >■ tib-l = ns — riL where n^ 
(or ub , depending on the interpretation of b as baryons or 
leptons, respectively) vanishes in the toy model. Within 



the standard Boltzmann approach, an evolution equation 
for ub-l can be obtained by subtracting the Boltzmann 
equations ^ for /;, and fj,. Dividing the left- and right- 
hand sides of Eq. ^ by p° = Ep and integrating over the 
momentum space one finds |32| : 



J_d_ 
a^dt 



{a'uB-L) = Jdnl [cr"(i + /,) - ct^h 

-cr\i + h) + ct-'h], (5) 



where a is the cosmic scale factor and t the proper time. 
Let us stress that, when inserting Eq. ^, the structure 
of the in-out matrix elements would lead to a generation 
of an asymmetry even in equilibrium. This inconsistency 
originates from a double counting of decay followed by 
inverse decay and scattering with a ipi in the intermedi- 
ate state. If the quantum statistical terms are neglected 
it can be removed by explicitly subtracting the on-shell 
part of the s-channel scattering, a procedure known as 
real intermediate state subtraction [3j, [SJ] . This is an 
example for double-counting mentioned in the introduc- 
tion, and can be completely resolved in the CTP ap- 
proach [H [11 H P • 

The processes which contribute to the right-hand side 
of Eq. ([5]) can be classified as source terms Sq , which ac- 
count for the generation of an asymmetry, or washout 
terms Wo, which tend to deplete the asymmetry. The 
subscript should remind the reader that, in the standard 
approach, only zero-order gradient contributions are in- 
cluded. In the hierarchical limit, M2 ^ Mi = M, the 
integrated Boltzmann equation for the B ~ L number 
density in the comoving volume can be cast into the 
form [35I-I37II 



dY, 



B-L 



dz 



So{z)~yVo{z)Yi 



B-L ■ 



(6) 



where Y — n/ s is the yield and z = M/T is the in- 
verse temperature normalized by the mass of the lightest 
toy- neutrino^. Performing the usual approximations, one 
finds from Eq. ^ that the source term is given by 



S^{z) 



■■KZ^D{Y^n-Y;i). 



(7) 



where n = T / H\t=m is the so-called washout parame- 
ter and the thermally averaged dilation factor is given 
by the ratio of two modified Bessel functions, 7d = 
Ki{z)/K2{z). Note that the structure is completely anal- 
ogous to the phenomenological case. 

Equations ([3]) also imply that the washout term is 
given by 



Wo(^) 



KZJB-L ■ 



(8) 



^ We have replaced the derivative with respect to the proper time 
t in Eq. (|6]l by the derivative with respect to the dimonsionloss 



inverse temperature using the relation -jv 
the FRW universe. 



valid in 



where jb-l = z'^Ki{z). Note that in a symmetric con- 
figuration, i.e. for Yb-l = 0, the contribution of the 
washout term in Eq. ([6]) vanishes. We will use this prop- 
erty to calculate the source term in the closed-time-path 
formalism later on. The yield Y^-^ (z) obeys an equation 
similar to Eq. ([5]), dY^^/dz ~ —kzjd (^i/^i — ^0')- 

In the strong washout regime, i.e. for large k, there is 
a well known asymptotic solution of this system. Up to 
an overall numerical factor 0(1) it reads 



770 = YB-L{t -^ 00) Oi 



Zf K 



(9) 



where Zf is the so-called freeze-out inverse temperature 
which is determined by the solution of nzf^B-L = 1- 



B. Closed-time-path approach 

Quantum corrections to the semi-classical Boltzmann 
approach can be studied using nonequilibrium quantum 
field theory techniques which rely on the closed time path 
formalism. In the remainder of this section we briefly 
review the closed time path approach, its relation to the 
Boltzmann approach, and discuss the gradient expansion. 

In general the B—L asymmetry ([¥]) is given by the zero- 
component of the expectation value of the corresponding 
quantum mechanical current operator: 



j^{x) = 2z {[V^b{x)] b{x) ~ h{x) V^h{x)) 



(10) 



Note that (•) = Tr(p-), where the density matrix 
p characterizes the system at some initial time iinit- 
The time-evolution of such in-in expectation values can 
be described within the closed-time-path or Schwinger- 
Keldysh approach. A useful quantity are the Wightman 
propagators, 

D> (x, v) = 5< (y, x) = {b[x)b{v)) , (11a) 

D<(x,y) = Dy{y,x) = {b{y)b{x)) . (lib) 

In terms of these, the B — L current can be expressed as 

j^,{x)^2iV.,, [D>{x,y)-D<{x,v)\ \y=^ . (12) 

The time-evolution of the Wightmann propagators is de- 
scribed by so-called Kadanoff-Baym equations, which are 
self-consistent Schwinger-Dyson equations formulated on 
the closed time path. In the limit of coinciding arguments 
they read [ill, [ll 

[n, + m^]D^{x,y)\y^^ = i j &'^z[J:<{x,z)D>{z,x) 

tinit 

-I]>(x,z)i?<(z,x)]. (13) 

The self-energies S]< and I]> in Eq. p^ can be inter- 
preted as generalizations of gain and loss terms, respec- 
tively [HJ, [i3l- At one- loop level, see Fig. [21 they read 




FIG. 2: One-loop contributions to the self-energies of the toy- 
neutrinos (a) and toy-baryons (b), respectively. 



9i9*jG'i{y,x)D^{y,x) 



S^(x,2/) .,.,^_,^,„,„_,_^ 



(14) 



In addition, it is necessary to consider corresponding 
equations for anti-particles, differing by D — >■ £) and 
S — > E in Eq. P^ . where the latter arc defined analo- 
gous to Eq. (HI]) . Let us stress that D and G*^ in Eq. P^ 
are the full non-perturbative Wightmann functions de- 
termined by the Schwinger-Dyson equations for the toy- 
baryons and toy-neutrinos, respectively. 

Using Eq. p^ we can next derive an equation for the 
divergence of the baryon current. It reads 

V^j^ix) = 2z (□, + m^) [Dy{x,y) ~ D^ix,y)] \^^^ 
= - / &'^z{j:<{x,z)Dy{z,x)-J:>{x,z)D<^{z,x) 

tinit 

- S<(x, z)Dy{z, x) + S>(a;, z)5<(z, x)] . (15) 

Formally, this equation describes the full quantum time- 
evolution of the expectation value of the B — L cur- 
rent starting from arbitrary (Gaussian) initial states at 
t = iinit- Therefore, it is a suitable starting point to 
derive non-equilibrium quantum corrections to the stan- 
dard results. Note that in the comoving frame in a ho- 
mogeneous FRW space-time j^^ = (?t-b-l, 0), i.e. 






B-L 



K dz 



(16) 



where K ^ — sdz/dt = szH = sH\t=m / z. Thus, 
Eq. (jlSp indeed constitutes a quantum generalization of 
Eqs. © and ©. 

Typically, there exists a separation of fast microscopic 
time-scales, e.g. fmic ~ ^/M, and the macroscopic evolu- 
tion characterized by decay and expansion rates, imac ~ 
1/r, 1/H. The crucial observation is that the former 
sets the scale for the variation of the two-point functions 
D{x, y) with respect to the relative coordinate, given by^ 
s = X — y , while the latter set the scale for the vari- 
ation with respect to the central coordinate, given by 
X = {x + y)/2. If there is a strong hierarchy between 
microscopic and macroscopic time-scales, which is typi- 
cally the case for thermal leptogenesis, it is possible to 
expand Eq. (jisp in gradients with respect to X, and keep 



only terms up to a certain order (see e.g. [3^). For that 
purpose, it is natural to express the correlation functions 
in the Wigner representation, which describes the "fast" 
variations along the relative coordinate s in momentum 
space, 

D^{X,p) = Jd'^se'P' D^{X + s/2,X- s/2) . (17) 

Then, the gradient expansion of Eq. (jlSp in the limit 
Unit — >■ — oo can be obtained using the general rela- 
tion [H 

J&'^zA{x,z)B{z,x) 

= Jdnpe-^0{AiX,p)}{B{X,p)}, (18) 

where dllp = d'^p/ {2Tr)*. The derivative operator 

<>{■}{■} = ^(4'^af - 9«9i')){.}{.} (19) 

generates an expansion in gradients with respect to the 
"slow" coordinate dxi^ (the superscript refers to the first 
and second argument). In particular, Eq. (|15p can be 
formally expanded in powers of <C>: 



25'^Jm = 2?^JmIo + 2?^jJi 



(20) 



As mentioned above, the standard Boltzmann limit 
is based on the zeroth order in the gradient expansion, 
which reads (see Appendix IXj) 



-D^J^lo = -JdUp e(po){sO<(l + h) 



E? 



> fb 

-S°<(l + /5) + S0/^}i?,. (21) 

Note the close similarity of Ec^s. (PT|) and ([S])- The first 
two terms on the right-hand side can be interpreted as 
gain- and loss terms for particles, respectively, while the 
last two terms represent gain and loss terms for anti- 
particles. Note that anti-baryons have negative baryon- 
number, hence the relative minus sign of the second 
compared to the first line. Within the CTP approach, 
the rates of gain and loss processes are described self- 
consistcntly by the self-energies I]> and f]> . Their struc- 
ture is fixed by the CTP formalism, which resolves the 
ambiguities of the standard Boltzmann approach men- 
tioned before. Note that in Eq. ((2T|) . we already intro- 
duced the distribution functions for baryons, /;,, accord- 
ing to the so-called Kadanoff-Baym ansatz. 



D<=hDp, D> = {l + h)D^ 



p ' 



(22) 



together with an analogous relation for anti-baryons in- 
volving /;,. Here Dp = Dy — £)< is the spectral func- 
tion, which has a Breit-Wigner-like shape. In the quasi- 
particle limit it reads 



See 1291 . 13811 for a proper generalization to curved space-time. 



Dp{X,p) = 27rsign(po) 6 (g'^'p^p, - vr?) 



(23) 



The corresponding self-energies in Wigner representa- 
tion, see Eq. ((ni), read [ll| 

S^(X,p) = -JdUkdU,{2T:)^S{k -q-p) 

xg:g,Gl{X,k)D^{X,q), (24a) 

t^iX,p) = -JdIl^dU,{27r)^Sik-q-p) 

xg,g*Gl[X,k)D^{X,q), (24b) 

where dWp = d^p/(27r)^, and G> and D^ denote the 
Wightman propagators of the toy-neutrinos and the 
baryons, respectively (see Appendix IBJ for more details). 
The Wightman propagator of the toy-neutrinos is a 
two-by-two matrix with non-zero off-diagonal elements. 
The latter describe mixing of the heavy fields and, in 
particular, carry information on the CP-violation in the 
system. For a hierarchical mass spectrum of the toy- 
neutrinos one has 11211: 



gvac ^ ^medQ^rp^ agrccs with the one obtained from the 
quantum-corrected Boltzmann equations derived in [111 . 
Il2l | and incorporates medium corrections Il3l . In the 
limit of an hierarchical mass spectrum jlll Il2| : 



e{k,T) 



^{^ + lBMEi) + ME2)] 



(29) 



where the second factor accounts for the medium effects 
and E^2) = ^[{M^ + |kp)^ ± (1 - 4m^/M^)^k\ cos0]. 
Note also that the structure of the CTP source term 
automatically guarantees that no asymmetry is produced 
in equilibrium. This means that it is inherently free of 
the double-counting problem and no explicit real inter- 
mediate state subtraction is required. Finally, when ne- 
glecting quantum statistical terms and medium correc- 
tions, assuming kinetic equilibrium for -01, and inserting 
the quasi-particle approximation (|23p . one recovers the 
standard source term Sq given in Eq. ([7|) . 



G| = £.g| + £;^|, 



(25) 



so that Eq. (|24|) can be recast in terms of the CP- violating 
parameters, 



-2lm{gj/gi)linei 



(26) 



and the diagonal components G> of the Wightmann func- 
tion. Using the quasiparticle approximation and inte- 
grating Eq. (|24p over the time-component of the four- 
momenta we obtain gain and loss terms similar to those 
in Eq. Q , see below. 

Let us now discuss the CP-violating source term within 
the CTP approach. As has been mentioned above the 
contribution of the washout term vanishes in the sym- 
metric configuration, i.e. for /& = /j, and only the source 
term contributes to the divergence of the baryon current. 
Comparing the integrated Boltzmann equation 1^ for the 
yield and Eq. (|16p suggests to define the CP-violating 
source term as 



S{z) = K ■ V'j, 



(27) 



where the subscript means that one evaluates the right- 
hand side of Eq. ([15]) for a symmetric configuration. Us- 
ing the zero-order gradient contribution (|2ip and insert- 
ing the explicit expressions for the self-energies ([M|) . we 
find for the source term in the CTP approach, to zeroth 
order in the gradients: 

So = K ■ P'^j^lo, sym 

^2K\gi\^JdUpdUkdUgeipo)i2n)''S{k-p-q) 
X gl\k)Dp{p)Dpiq) ■ e{k,T) [U,{k) - 4^(fc)] 

X ([1 + /rb)][i + n'm - n'{p)n'{q)) , (28) 

where f,p-^ and QY^ are the distribution- and spec- 
tral function of the lightest toy-neutrino, respectively, 
and we have assumed that the baryons are in equilib- 
rium. The effective CP-violating parameter e(fc, T) = 



III. ADDITIONAL SOURCE TERM FROM 
HUBBLE EXPANSION 

In this section, we derive and discuss the additional 
CP-violating source term for the B — L asymmetry, which 
is generated due to the time-dependence of the tempera- 
ture of the thermal bath, f = -HT. 

An important feature of the CTP approach is that the 
CP-violating source term ([27| automatically vanishes in 
complete thermal equilibrium [111 Il2|. Thus, it must be 
possible to express it in terms of quantities which are non- 
zero only out of equilibrium. Obviously, the deviation of 
the toy-neutrino distribution from the equilibrium one, 
^fipi = fipi ~ J'^ 1 is such a quantity. As discussed above, 
a second one is given by T = dT j dt. Thus, one may 
decompose the source term ([27[) according to 

S{z) = K-V^j^Uym 

= 5(z)|a/+5(z)|j,+5(z)U^,j,. (30) 

The first contribution is given by Eq. (|28p , of zeroth order 
in the gradients. 



S{z)\Af^Siz)\o^K-V^j^\o, 



(31) 



The second contribution appears only at first order in the 
gradients. It can be computed by setting all species into 
local thermal equilibrium (LTE) (i.e. setting A/^j^ — 0) 
with time-dependent temperature T(i), 



Siz)\^ = Siz)\i^^^K-V^j,\i; 



sym 



(32) 



Finally, the third contribution requires both Af^-^ ^ 
and T 7^ 0. In the remaining part of this section, we will 
show that a source term proportional to the expansion 
rate, S oc H, indeed follows from S{z)\f. For a discussion 
oi S\^f^rp we refer to section HVl 



The first-order gradient contribution to the full ex- 
pression for the B — L current, Eq. P^ . can be ob- 
tained straightforwardly using Eq. (jlSp , and keeping lin- 
ear terms in the derivative operator ( "Poisson bracket" ) 
defined in Eq. ([TO| (see Appendix |X| , 

+2<}{i:%}{Dh} + 2<}{El}{Dp} 
~20{t%}{D^} - 20{f:U{DF}) . (33) 

The self-energies Si. are given by Eq. (|24| with the first- 
order gradient solutions of the Kadanofl-Baym equations 
for the toy-neutrino propagator, 

G'^ = -^[o{gR,ff^,gA} 

where we have introduced ^{A, B,C} = <;^{A}{BC} + 
A<){B}{C} (see Appendix [B]). Furthermore, Dp = 
(!)> -|- -D<)/2, with analogous relations for self-energies, 
and Dh = {Dii+DA)/2, involving retarded and advanced 
functions (see Appendix [XI for more details). 

As discussed above, the source term S{z)\f can be ob- 
tained by evaluating Eq. (|33p for a symmetric system in 
local thermal equilibrium. After a somewhat tedious cal- 
culation, for which we refer to Appendix O one finds 
that at leading order in the toy-neutrino coupling, and 
in the hierarchical limit, it is given by 



S{z[ 



I ■ = K ■Vf'i 1^^^ 

\T ~ JmII, syrn 

= 2K Igil^JdUpdllkdUg e{po){27r)^S{k - p - q) 
X gl\k)Dpip)Dpiq) 
xe^ik,T){f;i{k)[l + r,''ip)][l + r,\q)] 

+ [i+/;!(fc)]/r(p)/r(9)}, (35) 



with 



<k,T) 



,T_ 
2T 



^ + ^u-d,\L,{k;T), (36) 



and 



Lhik;T) =167r/dn, 6(4) Dp{i) [i + /;;"{£)] 
-V -V 



(fc-£)2-m2 (fc + £)2-TO2 



, (37) 



where u ~ (1,0,0,0) in the co-moving frame, /^"^(fc) = 
l/[cxp(ii-fc/r(t)) — 1], and P denotes the principal value. 
This additional source term is the main result of this 
work. Note that it is proportional to the expansion rate, 
since T = —HT, as expected from a gradient contri- 
bution. Its structure is qualitatively different from the 
zero-order source term (|28[) , since it can be non-zero even 
when the heavy toy-neutrinos are in equilibrium with the 
thermal bath. We stress that the structure of this source 



term is, nevertheless, in agreement with the Sakharov 
conditions: it is proportional to e'""'^, i.e. it requires CP 
violation, and it requires a deviation from equilibrium, 
which is, however, described by T here. 

Let us now discuss the implications of this source term. 
An explicit expression for ef can be found in Eq. (|C24p . 
In the strong washout regime most of the asymmetry 
is generated at temperatures T < M. In this limit we 
obtain for the new CP-violating parameter 



j,(fc,T)^-6™'= 



H 

2^ 



(38) 



Inserting H 
we find 



IMy^T^/Mpi with .g* = 106.75 [3C 



.{k,T) 



-2.7 e^ 



T 



M, 



PI 



(39) 



This is clearly much smaller than the CP-violating pa- 
rameter t""''^ itself for typical temperatures T ~ lO^GeV. 
However, one should keep in mind that the structure of 
the source term (|35p differs from the usual one. In par- 
ticular, if the neutrino is very close to equilibrium, the 
zero-order source term (P5|) is suppressed, which could 
partly compensate the smallness of tf. This is typically 
relevant in the strong washout regime. Therefore, we 
expect that it is legitimate to neglect the gradient cor- 
rection psp . unless for extremely strong washout. In or- 
der to estimate its impact, we apply the approximations 
common for the standard approach to Eq. (|35p . and use 
Eq. ([55]) . In this limit Eq. ([55]) can be re- written as an 
additional source term to Eq. ([5]) , 



^^^ = 5o(z) -t- Sf{z) - Wo{z)Yb-l 



where Sf is given by 



Sf{z) 



-K-lB-L- 



n \t=m £ 

M 27r3 



(40) 



(41) 



Note that we neglected gradient corrections to the 
washout term here. This can be justified by the obser- 
vation that the latter have the same qualitative struc- 
ture than the zero-order washout terms, unlike the source 
terms. 

Let us assume now, that k is very large and the sys- 
tem is very close to thermal equilibrium. In this limit 
Y^ ~ y^' and the standard source term can be neglected. 
Thus, the rate equation for the baryon number density 
in the comoving volume takes the form: 



dY, 



B-L 



dz 



Sf{z)~Woiz)YB-i 



(42) 



Using the method of steepest descent we obtain an ap- 
proximate solution for the asymptotic value of the gen- 
erated asymmetry: 



Vgrad = Ys-lit ^ OO) (X 



M 



Zf Mpi 



(43) 
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FIG. 3: Ratio of the final asymmetries rj = Ys-Lit — > oo) 
obtained from the numerical solutions of the rate equations 
with both the gradient and the zero-order source terms taken 
into account to that with only the zero-order source term 
present for various values of the washout parameter k and the 
toy-neutrino mass M. In the limit k — >■ oo the dependence is 
\riarad\/r)o oc kM/Mpi. 



where the freeze-out temperature Zf is determined by the 
same equation as in the zero-order calculation. Compar- 
ing this result with Eq. ^ we see that the asymmetry 
is not suppressed by the washout factor k but on the 
other hand it is strongly suppressed by the ratio of the 
right-handed toy-neutrino mass to the Planck scale. The 
numerical analysis confirms that for large values of the 
washout parameter Eq. (|43p indeed well approximates the 
exact result. 

In Fig.[3]wc present the ratio of numerical solutions of 
Eq. (j40[) to those of Eq. ^ . For a wide range of k and 
M the gradient terms are subdominant and can safely be 
neglected. However, in the ultra-strong washout regime 
and for very heavy Majorana neutrinos they become the 
dominant source of the asymmetry. 

Important is not only the relative size of the two con- 
tributions but also the absolute value of the generated 
asymmetry. As is evident from Fig. 21 although for large 
K the contribution of the gradient source term dominates, 
the efficiency of leptogenesis in the standard cosmologi- 
cal setting is too small to reproduce the observed value 
of the asymmetry. However, the relative importance of 
the gradient terms could be strongly enhanced if other 
non-equilibrium phenomena would occur simultaneously 
with leptogenesis. Essentially, the size of gradient terms 
is determined by the scale of temporal or spatial inho- 
mogeneities. Since leptogenesis is often regarded to occur 
shortly after reheating, it is possible that non-equilibrium 
fluctuations occurring for example during (p)reheating 
are still present. Another possibility would be that a 



FIG. 4: Ratio of the generated asymmetry to the CP-violating 
parameter obtained from the numerical solutions of the rate 
equations with both the gradient and the zero-order source 
terms taken into account. In the upper right part of the plot, 
the gradient contributions dominate. 



phase transition occurs at the temperature relevant for 
leptogenesis, possibly related to the breaking oi B — L 
symmetry. If the phase transition is of first order gra- 
dient terms could have a great impact, similarly as for 
electroweak baryogenesis. Finally, we note that the anal- 
ysis presented here also applies to alternative mechanisms 
such as GUT-scale baryogenesis. Due to the much higher 
mass scale of ~ lO^^GeV, the gradient corrections are 
also larger in this case. 



IV. GRADIENT CORRECTION TO 
CP-VIOLATING PARAMETER 

In this section wc discuss the last term in the expan- 
sion (j30p of the source term. As we will sec, it can be 
interpreted as a correction to the CP-violating parame- 
ters ei. As has been argued above, ct are generated by 
the off-diagonal components of the Wightmann propa- 
gators of the toy-neutrinos provided that the latter are 
complex-valued. The first-order solution (p4| does have 
an imaginary component. Thus, it also contributes to 
the CP-violation in the system. Particularly interesting 
arc the following two terms of Eq. (|M)) : 



-^[G^O{nA}{GA} + 0{^fl}{^fl}^g 



(44) 



because they satisfy the condition (|25)) . To evaluate them 
we need explicit expressions for the self-energies and the 
propagators. In a state with zero (or small) asymmetry 
the self-energies corresponding to Fig. [2^ read |12| 






^t 



± -H'-' 



^[g*gj+9.g*)[Lh + hLil (45) 



where the functions Lh and Lp are defined in Ap- 
pendix |B] It has also been shown in [I4I that the di- 
agonal retarded and advanced propagators may be split 
into two real- valued diagonal matrices: 






(46) 



where the ofF-shell diagonal propa gato r and the spectral 
function arc given respectively by [IJ, [23] 



0^^ = - 



M? 



nj/ 



(fc2-Af2-n^^^)2 + i(n^pj) 



33 \2 



Q30 ^ 

^P IU2 



n« 



-^ nj/)2 + i(n^,^)2 



(fc2 - ap 
Substituting Eqs. (gS]) and (gHl) into Eq. (gl]) we find 



(47a) 
(47b) 



im£, = -0{nn{g//} 



In the hierarchical case to leading order: 






(48a) 

(49a) 
(49b) 



The decomposition coefficient e^ must be evaluated on 
the mass shell of the corresponding quasiparticle. Since 
we assume strong hierarchy of the masses here, Q^^^ 
evaluated on the mass shell of the i'th quasiparticle is 
negligibly small for i ^ j. Thus, the contribution of 
Eq. (|49bp can be neglected. Evaluating the contribution 
of Eq. (|49ap we obtain for the gradient correction to the 
CP- violating parameter:'^ 



^AfxT 



2[fc°P„L^] 



Mf 



(50) 



In the strong washout regime the toy-baryons are very 
close to equilibrium so that Lh{X,p) depends on time 
only through the dependence of the temperature T on 
time. Using once again the relation T = —HT we find 



k^'V^Lh ^ {E§-^ - Hk"^) Lh 



HT[{ 



E 



d 



T2 0(1/T) T dE 



\Lh] =HT-F, (51) 



where E and k are the energy and momentum of the de- 
caying toy-neutrino and P is a dimensionless function of 
T, k and Mi. In agreement with our expectations the 
correction to the CF- violating parameter is proportional 
to the Hubble parameter H. Evaluating Eq. ([SH]) in the 
radiation dominated universe we obtain in the hierarchi- 
cal case : 



Ml M^ 



^AfxT 



M, 



PI 



M| 



73 



(52) 



^ This effect exists also for the vertex contribution to the CP- 
violating parameter. However, there it is of a higher order in the 
couphng constants and for this reason we do not consider it here. 
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FIG. 5: Dependence of the function T on the dimensionless 
inverse temperature for the momenta of the decaying particle 
in the range 0.25 < |k|/T < 4. 



where we have introduced F = 2 ■ 1.66(7* ■ F for conve- 
nience. Let us analyze Eq. (|52p term by term. Since most 
of the asymmetry is generated at T ~ Af 1 , at the epoch 
of leptogenesis z ^ 1. The dependence of the function 
T on the dimensionless inverse temperature z for vari- 
ous momenta of the lightest toy-neutrino is presented in 
Fig. [5] Roughly speaking JF ^ 0(10), i.e. it is of the 
order of 2 • 1.66^^. For hierarchical mass spectrum, 
M1/M2 ^0.1, the requirement of successful leptogenesis 
implies that the mass of the lightest right-handed neu- 
trino should lie in the range Mi = 10*^ — 10^^ GeV. Even if 
we take the larger value, Afi ^ 10^^ GeV, its relative size 
compared to the Planck mass, Mpi = 1.2 • 10^^ GeV, is 
still very small, of order of 10~^. Consequently, the new 
contribution to the CP-violating parameter is strongly 
suppressed, primarily due to the smallness of the mass 
of the lightest right-handed neutrino as compared to the 
Planck mass. Let us note again that in GUT baryoge- 
nesis scenarios where the generation of the asymmetry 
takes place at higher temperatures, T ^ Mqut, the rel- 
ative suppression of the new contributions would be less 
pronounced. 



V. SUMMARY AND OUTLOOK 

In this paper we have derived gradient corrections to 
the kinetic equations for leptogenesis, that are neglected 
in the standard Boltzmann approach. 

We have found that there is an additional CP-violating 
source term Sf with a qualitatively new structure, which 
arises due to the time-dependence of the effective tem- 
perature of the thermal bath. It does not vanish even 
if all particle species are in local thermal equilibrium. 



For a standard cosmological background it is compara- 
ble to the conventional one only if the washout param- 
eter K ^ Mpi/Mi and can therefore be safely neglected 
in standard thermal leptogenesis. However, it becomes 
dominant in the limit of ultra-strong washout and for 
very heavy Majorana neutrinos, and can play an impor- 
tant role for alternative baryogenesis mechanisms oper- 
ating at very high scales, e.g. at the GUT scale. 

We have also analyzed a contribution to the effective 
in-medium CP-violating parameter which is induced by 
the gradient terms. Just like the Sj, term, the new con- 
tribution to the CP-violating parameter is suppressed by 
the small ratio of the heavy neutrino mass to the Planck 
scale. 

The gradient terms could be greatly enhanced in a 
non-thermal environment, where large temporal or spa- 
tial gradients can occur. This may be relevant when the 
reheating temperature is very close to the right-handed 
neutrino mass, as is often required to avoid the overpro- 
duction of gravitinos in supersymmetric scenarios. An- 
other conceivable situation is that the seesaw scale is as- 
sociated to the breaking of a symmetry, possibly B — L, 
in which case a phase transition could occur at tempera- 
tures relevant for leptogenesis. If it is of first order, gra- 
dient terms can play a major role, similar to electroweak 
baryogenesis. The additional gradient source term could 
then even allow to lower the scale of leptogenesis without 
having to rely on resonance effects. This is left for future 
work. 



Using analogous definitions for D and the self-energies, 
Eq. (jlSp can equivalently be written as 

P^j^x) = ij&^z e(xO - zO)e(z" - t„„t) (A2) 

x[l]i?(a;, z)Dp{z, x) — Sp(a;, z)Df{z, x) 

~I]f{x, z)Dp{z, x) + f^p{x, z)Dp{z, x)] . 

Note that we have expressed the integration limits in 
terms of the usual 0-function. It is helpful to absorb 
Q{x'^ — z*^) into retarded and advanced propagators, 

DR{x,y) = e(x"-y")7^p(.T,y), 



DAix,y) 



-e(2/°-.T")i?p(x,2/), 



(A3) 
(A4) 



again with analogous definitions for D and the self- 
energies. The Wigner representation of the various two- 
point functions reads 

A{X, p) = {-iY^d'^s e'-P' A{X + s/2, X - s/2) , (A5) 

where p = 1 for A = D^, Z)p, Tip^Tip, and zero otherwise. 
Using Eq. ([T51) and assuming tinu — > — oo. the Wigner 
representation of Eq. (|A2[) is 

'D^J^.{.x) = - ij'dUp (A6) 

X [er'0{T.F}{DA} + e-'^{E^}{I?^} 



-iO 



{tF}{DA}-e-'<'{tR}{DF}\. 



All two-point functions are evaluated at the point (AT, p) 
in phase-space. The gradient expansion formally follows 
by expanding the exponentials in powers of the derivative 
operator ex i9jf defined in Eq. (fT9)) . 
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Appendix A: Gradient Expansion 

In this appendix we derive the gradient expansion of 
the CTP evolution equation Eq. (fT^ for the B — L cur- 
rent, following the lines of |22l. l39l - l4ll |^. For that purpose, 
it is convenient to switch from the Wightman functions 
to an equivalent representation in terms of the so-called 
statistical propagator Dp and the spectral function Dp, 

Dp{x,y) = \[D>{x,y)+D^{x,y)], (Ala) 
Dp{x,y) = i[Dy{x,y)-D^{x,y)]. (Alb) 



* For simplicity, we will work in fiat space-time here. The relevant 
equations can be easily generalized to FRW space-time using the 
results of ref. 1291 . 



1. Zeroth order 

Evaluating Eq. (|A6p at zeroth order in the gradient ex- 
pansion yields 

-tPpDA - f^lDp} . (A7) 

By substituting p — > —p in the second line and using 

D^F^A,B.^^){X,p) = ^(F,fl,A,5)(X,-p), (A8) 

Dp{X,p) = -DpiX,~p) 

= -i[DniX,p)-DA{X,p)], 

this equation can be simplified to 

'D^J^.{x)\o = -/ cUlp [j:%Dp - J:^pDf] . (A9) 

Using Eq. (jAlap to re-express this result in terms of 
Wightman functions (note factors i in Eq. (jASp ) yields 

vr^^.{x)\o - -Jdnp [j:%d> - ^Id<] ■ (Aio) 

Finally, inserting the identity 1 = 0{po) + 0(— po) a-nd 
using Eq. (jASp yields, 

-S2,5> + E^5<] , (All) 
which corresponds to Eq. pT|) . 



10 



2. First order 

The first-order gradient contribution to Eq. (|A6[) con- 
sists of two parts. One obviously involves the linear 
term in the expansion of e~'^. In addition, it is impor- 
tant to realize that the self-energies, see Eq. p^ . con- 
tain an internal toy-neutrino line described by the out- 
of-equilibrium propagator G*^ . In order to obtain a con- 
sistent gradient expansion, it is important to perform also 
a gradient expansion of the equation of motion for G*^ , 
which is discussed in Appendix |B] This means we also 
have to expand the self-energies. 



E = S'^ 



Si 



(A12) 



Doing similar manipulations as for the zeroth-order 
equation yields 



V^^j 



/ill 



'/ (Mp \j:\Dy - T.lD< (A13) 



where Dh{X,p) = [Dr{X,p) + DAiX,p)]/2. Inserting 
again 1 = 9(po) + ©(^Po) and using Eq. (|A8|) yields 



Appendix B: Gradient expansion for real scalar 
fields 



The dynamics of the system of real scalar fields is 
described by the non-equilibrium generalization of the 
Schwingcr-Dyson equation [12l | 



{G-'y^{x,y) = {^-'yH^,y)-^''{^,y): 



(Bl) 



where G'-' is the full dressed propagator of the "heavy 
neutrinos" , ^*-' is the diagonal propagator of the free 
fields and 11'-' is the self-energy. Let us now split the self- 
energy matrix II into the diagonal, II , and off-diagonal, 
n , components and introduce a diagonal propagator Q 
defined by the equation 

g-'ix,y)=^-\x,y)-nix,y). (B2) 

Subtracting Eq. ((B2|) from Eq. (|BT|) we find: 

G-\x,y)^g-\x,y)-n{x,y). (B3) 



Multiplying Eq. ((B3)) by G from the left, by G from the 
right and integrating over the closcd-time-path contour 
[42 . |43| we obtain a formal solution for the full non- 
equilibrium propagator [I^l : 

G^ ix, y) = Q^ (^, y) - JJ^^u^^v e{u'>) eiv°) 

X [Gr{x, u)n^ {u, v)gA{v, y) 
+G-^{x, u)n^{u, v)gA{v, y) 
+GR{x,u)nR{u,v)g^{v,y)\. (B4) 



If the mass spectrum of the heavy scalars is strongly 
hierarchical, i.e. M^ <C M2, one can approximate the 
full propagators G on the right-hand side of Eq. (jB4p by 
the corresponding diagonal propagators Q. That is, in 
this approximation the dynamics of the diagonal and off- 
diagonal components of G is completely determined by 
the dynamics of Q. Wigner-transforming the resulting 
expression we obtain to leading order in the gradients 
Gi 



G>< 



G% 



r>, where 



G| ~ g^^ [GrJI'^Qa + Q^TIaCa + CRflRQ^] (B5a) 

G\c^-i[(>{gR,ff^,gA} 

+^{g^,nA^gA} + ^{QR,nR,Q^}\ ■ (B5b) 



The generalized derivative operators in Eq. (|B5bp arc de- 
fined by 



0{A, B, G} = (>{A}{BC} + A<}{B}{G} . 



(B6) 



A system in thermal equilibrium is stationary. There- 
fore, in equilibrium the right-hand side of Eq. (|B5b|) van- 
ishes. One can expect that in the early universe this 
contribution is proportional to the expansion rate H of 
the universe. A substitution of Eq. (jB5bp into the ex- 
pressions for the Wigner-transforms of the self-energy, 
Eq. ([M|) . gives us the first-order corrections S> to the 
self-energies. 



Appendix C: Derivation of the additional source 
term 



In this appendix we derive the additional GP- violating 
source term Sf{z) . Our starting point is the CTP evolu- 
tion equation for the B — L current P^ . The additional 
source term arises at first order in the gradient expan- 
sion, sec Eq. (jAlSp . Inserting baryon-symmetric propa- 



gators, D>,D- 



r^sym 



as well as the self-energy given 



in Eq. ^ into Eq. (|Al3)) yields 

S{z)\l=K■V^'j^{x)\l,sym^ 

2K lmig*gj)J dUp J dU^ J dUke{po){2TTfSik ~ p - q) 
X [lm{Gl;'')D'^"'D'^"'-Im{G]i'')Diy"'D'^"']. (CI) 

Here we have neglected the contributions in the second 
line of Eq. ()A13|) , similar as in [23 , and the first-order 
toy-neutrino propagators G^ are given by Eq. (jB5bp . For 
brevity, we suppress the superscript of D'^^^^ in the fol- 
lowing. 

Next, we set all particle species into local thermal equi- 
librium in order to obtain Sf = iSjf"^^. This means 
that the propagators fulfill the Kubo-Martin-Schwinger 
(KMS) relation 



D 



LTE 



(X,fc) = e^('=-")i?^^^(X,fc), 



(C2) 
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where u = (1, 0, 0, 0) in the co-moving frame of the ther- 
mal bath, and (3 = I3{t) = 1/T{t). Note that this condi- 
tion impUes that fb = l/[e^('=") - 1], see Eq. ^. For 
the toy-neutrino, the condition of local cquihbrium im- 
phes that the zero-order propagators G> fulfill KMS. In 
particular, this implies that 

g^^^(X,fc) = e^('=-")g^^^(X,fc). (C3) 

Again, in the following we suppress the superscript. 
These relations greatly simplify Eq. (|C1|) after inserting 
Eq. (jB5b[) . As an example, we consider the following con- 
tribution: 

= Pfc.g!!)Z)<D< - {dk^g^)DyD>] X v^iTig'/ 

-[v^{g^)D<D< - {v^g'^)DyDy] X dk^iv'ig'i 
+[g'^D^D^ - g^^DyDy] X <>{n'i}{g'i} . (C4) 

Using the LTE relations, we see that the third line can- 
cels. Furthermore, using again LTE, we find 



dk^gt^ = dk4e^^''^'5'^] = I3uc.e^^^-^'g 



f3{k-u)^ii 



+ e'3(*'-"'afc»^i3, (C5) 

_^ g/3(fc.n)pag«: ^ (C6) 

Inserting these relations in the first and second line of 
the above equation, we find that there are again some 
cancellations, 

0{gi:,ivig'2}D<D^ - 0{gii,ii'ig';i}DyDy = 

= g'^D^D^ {{f3u ■ V - (I?"/3)fc ■ud^.){IV'ig'i)] . (C7) 

Similar simplifications are obtained for the other two 
first-order gradient terms of G>*-',see Eg. (jB5bp . To 
shorten the notation we introduce the differential opera- 
tor 



pa^pa_ Ka^X)<Tp^g^^ 



(C8) 



Then, putting everything together, yields 

= -Klm{g*gj)JdnpdUkdUg {2nfd{k -p- q)Q{po) 
X Pu4lm[D'^{ilVlg'J)]g';D<D^ 
-Im['D'^{tg^lV^)]gi^D^D^ 

+ im{igl'b'-g'i - i'b''g'^g^i)i^iD^D^] . (C9) 

The parts where the space-time derivative acts on a re- 
tarded or advanced propagator gR[A) are suppressed, 
since they only depend on temperature via the thermal 
mass, which is of higher order in the coupling within the 
toy-model. Therefore, we neglect the third line. In ad- 
dition, using relations between retarded and advanced 



quantities and interchanging i <->■ j in the second line 
gives 

5(z)| J, = 2/v Ig.p/dnpdHfcdn, (27r)4<5(fc - p - g)e(po) 
x^f{k,T)[giDyDy+g';D<D4 , (CIO) 

where 

e^(fc,T) = \lm{gj/g,)liu^'b^lu,[ig^i{k)Ji^^{k)].{Cll) 

In the hierarchical limit, one obtains for the lightest toy- 
neutrino (i = 1): 

e^(fc,r) = ^f\k,T) = -e"'^- ^uj)^ LH{X,k), {G12) 

where e"""^ is the CP-violating parameter in vacuum. The 
loop integral Lh parameterizes the real part of the self- 
energy of the heavy real fields, see jl2 |. 

Lh{X, k) = 16^/dn, DF{q\ T)Dh{k - q) 

= Lh{k;T{t)) . (C13) 

In the co-moving frame: 

u^V- =T>'=V°- ^k'd,o - T A + |A:°a,o 



-H 



T^ + k"d,o 



(CM) 



This result can easily be generalized to a frame that is 
boosted with respect to the co-moving frame: 



u„2?° 



-H 



T^ + {u-k){u-dk) 



(C15) 



Thus: 



e^(fc,T) = e-x^ 



T^ + {u-k){u-dk) 



xLhik;T). (C16) 

In the quasi-particle limit, the loop integral is given by 



[\ + h{q-u-T) 



(fc -|- qY — m^ 



(C17) 



where q^ = m?,q^ = ^Jm? + q^ and fb~fi = f^^ for 
an approximately symmetric medium. The temperature- 
derivative can be evaluated using the result in [l2[: 



d 1 f°° 



X In 



T-^fnE;T) 

(2i; + |k|)2-fc2 



(2S-|k|)2-fc2 
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where fc^ = M^ and to « 0. The derivative with respect 
to the momentum is given by 



d 
—Lh[k;l) = iOTTjdi.^ 

2{k-q) 2{k + q) 



Lhik;T) = 167:JdIil [i + .^^^(g • u;T)] 



xV 



[{k — q)^ — rn^J^ [{k + q)^ — rn^ 



(C18) 



This expression is covariant, i.e. vaUd in any frame. In 
the co-moving frame, where u = (1,0,0,0), we need to 
compute only the derivative with respect to A:°. In a 
general frame, this corresponds to u ■ dk- Thus, 

d 



{u ■ k) X u ■ —Lh{k-T) 



k UkoLfi[k] 1 )\cQmoving 

2u ■ {k ~ q) 2u- {k + q) 



xV 



[(fc - g)2 - ?7l2]2 [(fc + g)2_j^2]2 



(C19) 



The above integral contains a "vacuum" and a "medium" 
part, where the latter is proportional to f^^{q-u; T). An 
explicit calculation shows that both integrals are well- 
defined (no UV, IR or on-shell-pole divergences), and 
yields (for m = 0) 



iu-k)iu-dk)Lr{k;T) 



In 



{u-k){u-dk)Lr''{k;T) 



Ek 
27r|k| 
Ek 



Ek - |k| 



Ek 



cBE 



* rif- 



Ar 



27r|k| 
AE^Ek 



|k|)^ 



M^-AE^{Ek + m)^ 



(C20) 



(C21) 



Here |k| is the momentum of the decaying particle in 
the comoving frame (rest-frame of the medium), and 
Ek = \/A/2 + k^ is its energy. For the part containing 
the temperature-derivative of L^, we use 

T-^f^iE; T) = -E^f^^{E; T) . (C22) 

Then one obtains 

T^ + {u-k){u-dk) Lh{k-T) 



Ek 
27r|k| 



In 



Ek 



Ek + |k| 



a/,^^ 



/^ OE 



In 



M" - AE^(Ek 



2E 
Ek 



In 



M4-4^2(^^ + |kp2 

(2^-K|k|)2-M2 



{2E - |k|)2 - M2 



(C23) 



It is easy to convince oneself that the remaining integral 
over the energy E is free of UV or IR divergences, and 
that the integral exists in the vicinity of the zeros of the 
arguments of the logarithm occurring inside the integra- 
tion region. Thus, we finally arrive at the result 



.(fc,r)-e"-x^ 



PH Ek 



27r|k| 



In 



Ek - |k| 



dE 



2E 



BE 



dE 

{2E^ 



In 



M^ 



Ek + |k| 
AE^iEk - |k|)2 



M^ - AE^Ek + IW 
'~M2 '' 



(2£:-|k|)2-M2 



(C24) 



In the non-relativistic limit |k| ^ T <C M, the main 
contribution to the integration over E comes from the 
region E <^ M , due to the exponential suppression in 
the Bose-Einstein function. Therefore, we may expand 
the logarithms in the above expression for |k|, i? <C M: 



In 


Ek |k| 
Ek + |k| 


^ -2|k|/ii;fe, 


(C25) 


{...} ^ lQE^\k\'/{M^Ek). 


(C26) 


Then it is easy to perform the energy integral: 




HdE^fE^^- 
Jo dE 


-J dEfr2E- ^^ . 


(C27) 



Using this, we obtain the leading contributions in the 
non-rclativistic limit. 



e^(A;,T)^e— X V — 



/3ff 
2 



1 

■K 



87rr2k2 
3i\/4 



(C28) 



The second contribution in the brackets is suppressed for 
T < M, so that we finally obtain 



<k,T) 



H 

2^ 



(C29) 
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